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A Monte Carlo simulation of the 0(4) At/)* theory in the broken phase is performed on a hypercubic lattice in 
search of an J = 1, J = 1 resonance. We investigate the region of the cutoff theory where the interaction is strong 
as it is there that a resonance would be expected to have a better chance to form. In that region the presence of 
an / = 1, J = 1 resonance with mass below the cutoff is excluded. 
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The simplest effective theory describing low en- 
ergy pion nucleon interactions is described by a 
chiral Langrangian, and involves the three pion 
fields TT"'^, TT^, TT^, the scalar field cr, the two nuclc- 
ons, their interactions, and the dimensionful pion 
decay constant /,r = 90 MeV. The a, although 
it has never really been seen, may exist in nature 
as a broad and quite heavy (« 700 MeV) isospin 
1=0, spin J=0 resonance. In that setting, it was 
natural to ask whether the presence of the p reso- 
nance was a consequence of the pion interactions 
or of some more fundamental interactions. 

Since the earliest such theory( Gell-Mann Levy 
sigma model ]l|) becomes an 0(4) Xcj/^ theory if 
the nucleon fields are neglected, this question re- 
duces to whether or not the 0(4) At/)'* theory can 
sustain an J = 1, J = 1 resonance in the broken 
phase. 

Today, in a very different setting, this same 
question is of interest again. The scalar sector 
of the Minimal Standard Model is also a four 
component A(/)* theory in the broken symmetry 
phase. The equivalent of the a resonance is the 
Higgs particle, the three pions are the Goldstone 
bosons, the pion decay constant is the weak scale 
fa = 246 GeF, and the scattering of longitudi- 
nally polarized vector bosons behaves exactly the 
same way as a scaled up version of tt-tt scattering 
(equivalence theorem 0). 

It is possible that the Higgs, like the a, is quite 
heavy and broad and may avoid detection at SSC. 
On the other hand, if the four component A(/>'* 
theory does indeed contain an / = 1, J = 1 
resonance, then since the p resonance, as it ap- 



pears in nature, is quite strong, its equivalent in 
the scalar sector of the Minimal Standard Model 
may have a good chance to be detected at SSC. 
Therefore, it becomes very important to know if 
another "signature" of the scalar sector, besides 
the Higgs resonance, is awaiting discovery at SSC. 

Earlier attempts Q to answer this question 
used analytical approximations such as Fade ap- 
proximants and found that an / = 1, J = 1 res- 
onance is present in the theory. It is not clear, 
however, whether or not the presence of this res- 
onance is an artifact of the approximations used 
and therefore the question of the existence of an 
/ = 1, J = 1 resonance in the theory has yet to 
receive a conclusive answer. 

We report here on a Monte Carlo simulation to 
shed new light on this old question [Q. The 0(4) 
A^"' theory has the lattice action 
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and was simulated on the lattice in the A — > oo 
limit. In that limit, the theory has the strongest 
interactions and a resonance probably has a bet- 
ter chance to form. 

The simulation was done on a hypercubic lat- 
tice A of spatial extension L and time extension 
Lt using an incomplete heat bath algorithm ||] 
on the CM-2 machine at SCRI. Since, on a fi- 
nite lattice, the direction of the symmetry break- 
ing changes from configuration to configuration, 
the same approach as in B was used to disentan- 
gle the massive scalar field a from the Goldstone 



modes. 

The mass rria of the a field was measured using 
standard techniques. To measure the energy of 
the lowest laying state in the / = 1, J = 1 chan- 
nel, an operator carrying these quantum numbers 
needs to be constructed. The simplest such oper- 
ator is: 
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where summation over repeated indices is as- 
sumed, a, 5, c are the isospin indices, tabc is the 
totally antisymmetric tensor, rh is the ra'th Eu- 
clidean unit vector of the time slice At , and m G 
[1,2,3] is the z-component spin index. Unfortu- 
nately, the time slice connected correlation func- 
tion of this operator gives a very weak signal. To 
get a better signal an operator that extends over 
several lattice spacings needs to be constructed 
using a trial wave function for the two-pion state. 
The "bag" and "bound state" meson wave func- 
tions were considered M . The latter gave a better 
signal with an affordable cost in computer time 
and was therefore used for the simulation. Us- 
ing as a trial wave function for the two pions at 
relative position i?, the hydrogen wave function 
^n=2,i=i,m(-R), an / = 1, J = 1 operator with 
total 3-momentum zero was constructed: 

Oc.™(t) = E E 



|i?|exp(-|i?|/2ao)n,™(e,0)7r>^^^ea6c 
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where i? is a three-dimensional cubic box centered 
at the origin and contained in At , ap is the "Bohr 
radius" in lattice units, 9 and </> are the azimuthal 
and polar angles of _R, and li.,„ is the I = 1 spher- 
ical harmonic (up to a multiplicative normaliza- 
tion constant). The parameter oq can be given 
any value. A very large oq will cause the expo- 
nential to decrease very slowly and then the sum 
over R will have to be carried out over a box B as 
large as At. Since this can increase the computer 
time significantly, a smaller oq has to be used 
so that the size of the box that contains the im- 
portant contribution from the exponential can be 
made smaller. However, oq cannot be made very 
small because the signal becomes weaker as ao de- 
creases. An optimal choice of aq and B was found 
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Table 1 

Energy E in the / = 1, J = 1 channel. 



to be oq = 2, and B extending from —3 to -1-3 in 
each of the three directions. The operator Oc,m{t) 
couples to the / = 1, J = 1 states. The energy 
of the lowest laying state in this channel can be 
found by looking at the time slice connected cor- 
relation function Cc,m(i) = < Oc,m(0)Oc,m(0* >c 
of this operator. 

The simulation was done on an L = 8, Lt = 16 
and L — 16, Lt = 16 lattice and for three values 
of the hopping parameter k, = 0.305, 0.310, 0.330. 
These values were chosen so that a wide range of 
TTio- will be covered (they also coincide with some 
of the values used in 0). The expectation val- 
ues of Ci,i(t) = \ Y.cm Oc,m{t), were measured. 
Ci,i(0 ^^s fitted for a few different ranges of t, 
and the resulting energies E, together with the 
X^ per degree of freedom for each fit are given in 
table 1 for the L = 16 and L = 8 lattices. 

The effective energy Eeff{t), for the two time 
slices at t-1 and it, is plotted versus t for the three 
values of k in figures 1 (L = 16) and 2 {L — 8). 
The values of t omitted from those plots had an 
Eeff{t) with very large error. 

In a two-pion 1=1, J = 1 state with zero to- 
tal 3-momentum, the lowest 3-momentum a pion 
can have has one component equal to ^ and two 
equal to 0. The next one has two components 
equal to ^ and one equal to 0. The energy spec- 
trum in the / = 1, J = 1 channel is therefore 
expected to contain levels with energies close to 
the energies of these states, but slightly differ- 
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Figure 1. Eeff{t) for the L = 16, Lt = 16 lattice. 



Figure 2. Eeff{t) for the L = 8, Lf = 16 lattice. 



ent because of the interaction. For the i = 16 
lattice, these levels have energies Eq ~ 0.78 and 
El ~ 1.09 respectively, and are denoted by the 
dotted lines in figure 1. From this figure, it is 
clear that the observed energy levels are very close 
to the levels of two free pions. In fact, for smaller 
t, the levels are close to Ei , and for larger t they 
are close to Eq. Because the free two-pion levels 
are not very well separated at L = 16, the ob- 
served levels are probably a mixture of the two 
lowest ones. In that sense, the energies in table 1 
for the L = 16 lattice are probably a mixture as 
well. The fact that the observed levels correspond 
to a two-pion state and not to a resonance is also 
greatly supported by the fact that these levels 
change only slightly from k ~ 0.305 to k = 0.330. 
After all, in that range ma- varies from 0.225 to 
0.91. Therefore, if a resonance is present it must 
have energy larger than « 0.78 and is either too 
heavy (for example, heavier than 1.09) to be ob- 
served with this method, or is "hiding" between 
0.78 and 1.09. That the latter is not the case can 
be seen by looking at the energy levels obtained 
for the L = 8 lattice. From figure 2 it is clear that 



there are no energy levels below « 1.4. In fact, 
since the two lowest free two-pion states are well 
separated in this lattice, the lowest energy level 
is clearly visible in this figure . It is true that the 
limited statistics give Epff{t) only up to t — 4, 
but because the lowest level is well separated from 
the next one, the correlated y^ fit gives a good es- 
timate of the energy of this level. The numbers 
given in table 1 are indeed very close to the lowest 
free two-pion energy Eq ~ 1.50 and do not seem 
to change much for the different values of k. 

From this analysis it is evident that if an / = 1 , 
J = 1 resonance exists for k. > 0.305, it is unphys- 
ical since it must have energy larger than the cut- 
off. Of course there is always the possibility that a 
resonance with energy below the cutoff does ex- 
ist, but, because the overlap of the operator in 
eq. 3 with an / =, J = 1 resonance may have 
been much smaller than its overlap with the free 
two pion states present in that channel, it was not 
possible to observe with the statistics and lattices 
used. However, it is unlikely that this is the case 
since the operator in eq. 3 is more appropriate 
for a resonance than for a free two particle state. 



For the theory of GeU-Mann and Levy of low 
energy pion processes the k > 0.305 region 
corresponds to the region where the a particle 
mass is greater than approximately 180 MeV (or 
equivalently the cutoff is less than approximately 
1.3 GeV). [| Since m^ is expected to be much 
larger than 180 MeV, it is clear from the results 
that the existence of the p resonance in nature 
cannot be accounted for by this low energy the- 
ory alone (in contrast to Q). 

For the scalar sector of the Minimal Standard 
Model the k > 0.305 region corresponds to where 
the Higgs mass is greater than approximately 
500 GeV (or equivalently the cutoff is less than 
approximately 3.5 TeV). The Higgs mass is of 
course not known but its upper bound is placed 
at around 650 GeV (hypercubic lattice triviality 
bound). Thus the existence of an / = 1, J — 1 
resonance can be excluded with confidence for 
values of the Higgs mass above « 500 GeV. For 
K < 0.305 (Higgs mass < 500 GeV) the strength 
of the interaction becomes weaker and hence it is 
safe to say that if a resonance could not form for 
K > 0.305, where the interaction is stronger, it is 
unlikely that it will in this region either. For this 
reson it was not deemed necessary to investigate 
the K < 0.305 region. A numerical simulation for 
K < 0.305 not only is not necessary, but it would 
also be very costly since larger lattices will have 
to be used (the correlation length becomes larger 
than approximately ^ ^25 — 4.5). 

It must be emphasized that these conclusions 
are valid only within the realm of the scalar sector 
of the Minimal Standard Model. It is of course 
still possible that the physics that enters at higher 
energies may be able to produce such a resonance 
in very much the same way the physical p particle 
owes its existence to QCD. This resonance if it 
exists due to some higher energy theory it would 
have energy determined by that theory. In fact, 
it is possible that the energy of this resonance 
is determining the cutoff energy of the Minimal 
Standard Model. 
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^ The connection with the physical scale /^ (or /q) is 
made through the renormalized coupling constant g^ = 
3m^/£B. The value of gn for this value of k was taken 
from 



